Frustrated phase separation in two-dimensional charged systems 
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We study phase-separation frustrated (FPS) by the long-range Coulomb interaction in two- 
dimensional electronic systems with emphasis to the case of a metallic and an insulating phase. 
In the mixed phase the system self-organizes in terms of mesoscopic inhomogeneities of one phase 
hosted by the other phase. We analyze the cases of circular drops and of alternating stripes. As 
a first approximation we consider the density inside each inhomogeneity as constant and in some 
cases we test the accuracy of this assumption by a more involved local density approximation. We 
find that the transition from the uniform phase to the frustrated phase-separated phase changes 
order depending upon its geometric arrangement. Contrary to what was found in three dimensional 
systems, there is no upper bound for the size of inhomogeneities. This difference stands on the 
different role of the long-range Coulomb interaction and screening in two and three dimensional 
systems. We conclude that two-dimensional systems are more prone to mesoscopic FPS. 

PACS numbers: 64.75.+g,71.10.Hf,71.10.Ca 



-a 
o 

> 

OS 

o> 

00 

o 
in 
o 

-I— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

A variety of strongly correlated charged sys- 
tems shows a strong tendency to phase separation 
(PS). 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 Often these system are 
two-dimensional (2D) like the two-dimensional electron 
gas (2DEG)2ii£ or quasi two-dimensional like cuprates 1 ^ 
or some manganites^i^iSii 7 ^ 1 ^ 9 " 

Advances in local probe techniques have revealed the 
mesoscopic nature of the phase coexistence. For example 
cuprates phase separate into superconducting like regions 
and poorly metallic regions on the scale of ~ bnm ~ 
10 lattice constantsi 1 ^* 1 ^*^ The regions have round like 
shapes indicating that surface energy has an important 
role. 

Mesoscopic phase coexistence has also been reported in 
manganites i 6 i 14 i 15 i 16 i 17 i 19 i 20 i 21 Here the situation is quite 
complex. Some systems have large insulating/metallic 
clusters with fractal like interfaces. 20 In this case phase 
separation is dominated by disorder effects stabilizing 
one or the other of the two phases very close in energy. 
On the other hand scanning tunneling spectroscopy in 
thin films have revealed filamentary and drop like metal- 
lic/insulating regions in the scale of tenths to thousand of 
nanometers with smooth surfaces indicating strong sur- 
face energies. In both cases percolation of the metal- 
lic regions is closely correlated to abrupt changes in 
transport .i^iSi Finally neutron and X-ray scattering has 
revealed much smaller clusters, on the nanometer scale, 
in the bulk, 6 ! 16 : 17 : 19 

Using a local probe in the 2DEG, Hani and 
collaborators^*^ have shown that close to a puzzling 
metal-insulator transitionSS* 2 ^ 2 ^ the system becomes 
inhomogeneous at a mesoscopic scale. In addition 
the compressibility close to the transition, departs 
sharply from the predictions of an homogeneous electron 
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gas^ta^^^^™ 

Mesoscopic inhomogeneities are generally expected in 
systems with a density driven first order phase-transition 
in the presence of long-range forces. This phenomenon is 
well known in a variety of systems 29 ranging from neutron 
stars^ to spinodal decomposition hampered by elastic 
forces 3 ^ and it is also related to the problem of domain 
formation in ferromagnetic systems 

In charged systems the phase coexistence phenomenon 
is hampered by the long-range Coulomb interaction 
(LRC). Indeed a macroscopic charge imbalance would im- 
ply an electrostatic energy cost which grows faster than 
the volume in the thermodynamic limit. Thus it has been 
proposed that the system should break in domains in or- 
der to guarantee large-scale neutralityi 4 i 5 i 7 i 8 i 30 i 34 i 35 i 36 i 37 
In this way the charge is segregated over some character- 
istic distance but the average density at large distances 
is constant. 

When the scale of the inhomogeneities is mesoscopic, 
one can perform a general analysis of FPS indepen- 
dently of the specific short-range interaction favoring 
PS, in the same spirit of the Maxwell construction 
(MC). 4 i 5 i 30 i 34 i 35 i 36 i 37 The specific short-range interaction, 
however, will account for different physical situations 
with different coexisting phases. Due to this mesoscopic 
hypothesis the inhomogeneities can be treated as charged 
classical objects. Their size and their relative distance are 
determined by the competition between the LRC inter- 
action and the interface boundary energy. These effects 
determine the total "mixing energy" i.e. the excess of en- 
ergy with respect to the unfrustrated PS. To characterize 
the degree of frustration one can define a dimensionless 
parameter A given by the ratio of the characteristic mix- 
ing energy and the characteristic energy gain due to PS. 

In 3D systems it was shown that the LRC interaction 
favors uniform phases £££^1 Indeed the coexistence re- 
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gion shrinks as A increases and the uniform phases are 
stabilized at densities where the ordinary MC analysis 
would predict a PS state. In addition the size of the in- 
homogeneities in 3D systems has been shown to satisfy a 
"maximum size theorem" that says that inhomogeneities 
can not have all linear dimensions much larger than the 
screening lengthi? 4 !?^? 1 ? 

In this work we consider mesoscopic frustrated phase- 
separation for 2D electronic systems, that is elec- 
trons confined to a plane but subject to the 3D 
Coulomb interaction in the presence of a rigid back- 
ground. Muratov, instead, has considered the case of 
a d-dimensional system immersed in the (i-dimensional 
Coulomb interaction^ As a first approximation we ne- 
glect disorder effects and concentrate on the thermody- 
namic behavior in the clean limit. We analyze in detail 
the coexistence between a compressible phase (named 
"B") and an incompressible phase (named "A") corre- 
sponding to the physically relevant case of phase sepa- 
ration between a metal and an insulator. The results 
are easily generalizable to two compressible phases^ A 
non-rigid background produces peculiar effects close to 
an energy level crossing of the electronic phases2^2ii£ 
and will be discussed elsewhere in the present context 4i 

We consider two different types of inhomogeneities: 
drops of one phase into the other phase and alternat- 
ing stripes of each phase. We assume a uniform density 
inside each inhomogeneity. This uniform density approx- 
imation (UDA) is relaxed in Sec. IIVI and tested in Ap- 
pendix [B] against a more involved local density approx- 
imation (LDA). UDA is shown to be quite accurate for 
the evolution of global thermodynamic quantities. 

As in 3D systems, the LRC interaction stabilizes the 
homogeneous metallic phase at densities where non frus- 
trated systems would show PS . 

In contrast to the 3D case, in 2D the density driven 
transition to the inhomogeneous state changes order 
depending upon the different geometric arrangement. 
Within our approximation, one has a second-order tran- 
sition to a droplet state while the transition to the stripes 
geometry is always first-order like. 

Moreover in 2D systems the size of the inhomogeneities 
is not limited by the screening length. As a result we find 
that 2D systems are more prone to mesoscopic FPS than 
the 3D systems. 



II. GENERAL ANALYSIS 

As in the 3D case we assume a separation of length 
scales. In many models on a lattice it is found that short- 
range forces tend to phase separate the electrons between 
electron rich and electron poor re gionsi^2iiM4ii5ii£iil 
We assume that the short range forces lead to a double- 
well form for the energy of an hypothetical uniform phase 
as a function of density. We call Ja and Jb the en- 
ergy densities close to each minimum, which define the 
bulk energy of the phases. Long-range forces prevent 



large scale phase separation leading instead to domain 
formation. As a simplified assumption we assume uni- 
form densities (nB,nA) for the inhomogeneities of each 
phase (UDA) with sharp interfaces (soft interfaces have 
been considered by Muratov^) . This charge distribution 
is compensated by a rigid background of density n. The 
UDA will be relaxed in section llVl 

The free energy per unit "volume" of the mixed-phase 
reads: 

/ = (1 - v) fA (tia) + v j B (n B ) + e e ; + e a (1) 

where v denotes the volume fraction of the B-phase 
(Vb/V), e e i represents the LRC interaction energy den- 
sity and e a is the "surface" energy density. (Here and 
below, in order to keep a common nomenclature with 
the 3D case, we call "volume" a quantity with units of 
length squared and "surface" a quantity with units of 
length). Due to charge neutrality the global charge den- 
sity of electrons has to compensate the charge density of 
the background en. This leads to the following neutrality 
constraint: 

n = (1 — v) riA + v Us 

To proceed further we assume specific geometries for 
the inhomogeneities. We will consider the competition 
among the two following geometries: i) drops (disks) of 
one phase surrounded by the other phase and ii) a peri- 
odic structure of alternating stripes of the two coexisting 
phases. 

For the drops we divide the system in cells of radius 
R c enclosing one domain of the B-phase with radius Rd- 
Similarly for the stripes the cell has width 2R C and con- 
tains a B-phase stripe of width 2Rd- The volume fraction 
is related to these characteristic lengths by: 

--(*)' 

where for 2D systems a = 1,2 for the stripe geometry 
and the drop geometry respectively. In the latter case 
the cells are slightly overlapping with: ttR^ Nd = V (Nd 
indicates the number of cells in the system while V is the 
total volume). For the stripe geometry V = 2R c LNd 
where L indicates the length of the stripes. The total 
surface energy per unit volume can be parameterized by 
a quantity a with dimensions of energy per unit surface 
(actually length in 2D). It reads: 

N d Zd 



where is the surface of the domain interface inside one 
cell with Yjd = 2it Rd and = 2L for the drops and the 
stripes respectively. For inhomogeneities of the B-phase 
hosted by the A-phase, e a can be written as: 

a a-l_ 

e* = — av « (2) 
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The operation: A ^ B, v ^ 1 — v was named "phase 
exchange" in Ref.|3 Within the UDA the energy should 
not change under this operation^ For the stripe geom- 
etry case the surface energy does not depend on v and 
Eq. J5J) already preserves the phase-exchange symmetry. 
On the contrary for the drop geometry Eq. (J2J) is appro- 
priate only for small v. In fact for intermediate volume 
fractions the drops should deform and for v ~ 1 the B- 
phase should represent the host in which drops of the A- 
phase are immersed. Eq. @, instead, violates this phase 
exchange symmetry because it involves drops of B-phase 
in both cases. We can define a "symmetrized" interpo- 
lating form that is correct at the two extremes [y ~ 
and v ~ 1): 



— a[v{\-v)\ 

tie 



(3) 



This is enough for our proposal because we can anticipate 
that drops are stable only in a narrow region close to 
!/~0 and z/~ 1 (c.f. Fig. HJ). 

The long-range Coulomb interaction energy is com- 
puted in the Appendix [S] by dividing the systems in 
neutral Wigner-Seitz cells. For stripes the electrostatic 
energy is computed numerically. 

The electrostatic energy contribution due to the inter- 
action between different cells can be neglected for volume 
fractions close to v ~ and v ~ 1. In these limits, in 
fact, the inhomogeneities are far from each other and the 
interaction becomes irrelevant. We show below that this 
is indeed a reasonable approximation in the full range 
of volume fractions. As discussed in Appendix El we 
expect this approximation to be even more accurate for 
the drop geometry. Therefore for further analytical com- 
putations we use the approximate expression derived in 
Appendix El for both geometries: 

e 2 8 s 

e ei = —{n B - n A ) 2 R c - \v (1 - i/)] 5 Drops (4) 
£0 3 



e et = —(n B -n A ) R c 2 [v (1 - v)\ [-\ogv(l-v)\ 
£0 

Stripes (5) 

Here Sq is the static dielectric constant. For the drop 
geometry the electrostatic energy has been symmetrized 
similarly to the surface energy. 

The bulk free energies f A and f B appearing in Eq. 
do not depend upon R c . The cell radius is therefore 
determined by minimizing the mixing energy e m = e e ; + 
e CT at fixed v. For the two geometries we have: 



<7E 



e 2 (n B - n A ) 



4 (1 - v) 
u{y) 



for the drops and 
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<7E 



e 2 (n B - n A ) 



u{v) 



(6) 



(7) 
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FIG. 1: Top: Approximate u functions parameterizing the 
mixing energy for the drops [Eq. ©] and the stripes [Eq. (1 101) 1 . 
For the latter geometry we also provide the numerically evalu- 
ated expressions taking into account the electrostatic interac- 
tion energy (see Appendb(XJ, and the LDA (see Appendix. |EJ 
in order to test our approximations. Bottom: Expanded scale 
for v ~ and v ~ 1 : the drop geometry introduces a lower 
mixing energy contribution. 



for the stripes. 

Notice that the exponent could be anticipated from 
dimensional analysis. For arbitrary dimensionality d of 
the system we have: 



Rc 



a e 



e 2 (n B - n A ) 



times a function of v. 

Once R c has been eliminated the mixing energy can be 
put in a common expression together with the 3D case 



where: 



a d V (n B - n A ) 
£0 



u(v) = -j=u(l-u) 



u(v) 



for 2D drops, while: 

u(v) = 2V2i/(l-v){-logu(l-iy)]i 



(8) 



(9) 



(10) 
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for 2D stripes. The corresponding functions in the 3D 
case can be found in Ref. yj. The function u {v) encodes 
all the information about the geometry and it is repre- 
sented in Fig. for the 2D case. 

We see from Fig. ^that only close to v ~ and v ~ 1 
the two geometries compete, since for intermediate vol- 
ume fractions the stripe geometry is stabilized. The in- 
terplay among the two states is therefore well described 
by the analytical approximation for the electrostatic en- 
ergy Eqs. Comparing the two curves for the 
stripe geometry we see that the use of the exact numeri- 
cal expression will only produce small changes in physical 
quantities for v ~ 1/2 but cannot change the qualitative 
behavior. Thus, our use of the analytical expressions 
(Eqs. (©,101) of e e i is justified. 

We minimize the total free energy per unit volume with 
respect to n B — tia an d v as m Ref ■ US and obtain the 
coexistence equations for both 2D and 3D systems: 



A*b - fJ>A 



Pb - pa 



n B 



n A 



u{v) 



- n A \ 2 ~* v(l-v) 



\n B 



[n + (n B - n A ) (1-2 v)] (ji B ~ Ma) + 



of phases with negative compressibility and its full dis- 
cussion will be presented elsewhere^i 

Since the insulator is incompressible and thus electro- 
statically inactive, we have to consider only the excess 
density of mobile electrons counting from the insulating 
state. Without loss of generality we set n A = and 
fA = for the insulator. In this way our density n has 
the meaning of density deviation from the pure insulat- 
ing phase. Also without loss of generality we consider 
?i > so that the insulating phase is at low density and 
the metallic phase is at high density but our results ap- 
ply equally well to the opposite case simply by changing 
n — > — n. With these conventions we are reduced to study 
the problem of the phase-separation between a metal at 
finite n B and "void" (playing the role of the insulator). 
Both phases are in the presence of a uniform background. 
In the present case the neutrality constraint reduce to: 
vn B — n, while minimizing the free energy in the phase 
separated state the two Eqs. (|llf) reduce to a single equa- 
tion for the B-phase pressure: 



Pb{ti b ) = p n 



(13) 



where the function p B (n B ) is a property of the bulk phase 
[c.f. Eq. jnjl] and 



r d-l „2 



e 2 (n B - n A ) \ d du (v) 
so / dv 



where 



(11) 



91 



A,B 



driA,B 

and the pressure is defined as: 

Pa,b = —Ja,b + MAS riA,B 



(12) 



We have neglected a possible dependence of a on density 
which contributes with additional terms (see Ref. |34j). 

Eq. I|ll(l generalizes the ordinary Maxwell construc- 
tion to the charged system case. The former ([ia = Ms 
and pa = Pb) is recovered when e — » 0. The first of 
Eq. Ijllfl states that the most dense phase has the lowest 
chemical potential. The difference of chemical potentials 
becomes particularly simple for 2D drops. In this case 
it is (within our approximations) constant and equal to 
±8a 1 / 2 e/(3e ) 1 / 2 . 



III. SEPARATION BETWEEN A 
COMPRESSIBLE AND AN INCOMPRESSIBLE 
PHASE 

Now we want to investigate an inhomogeneous state 
for a 2D system which involves an incompressible A- 
phase hereafter "the insulator" and a compressible B- 
phase hereafter "the metal" ~ 

This will be appropriate for example for a doped Mott 
insulator. The case of the 2DEG requires consideration 



Pn 



a e 



n B 



du (y ) u (i>) 



Eq. (|13fl with the neutrality constraint determines the 
behavior of the local density and volume fraction as a 
function of n. 

In the limit of e — > we obtain the equation pb^b) = 
which corresponds to Maxwell construction for the case 
of phase separation between a self-bound neutral fluid 
and vacuum. 

When e /: 0, p m represents the pressure due to the 
presence of the long-range Coulomb interaction and the 
surface energy, that we label " mixing pressure" and must 
be balanced by the B-phase pressure. 

For the drop geometry one has: 



Pn 



V3 



a e 



(14) 



which depends on the global density of the system only. 

Instead, in the stripes geometry the mixing pressure 
depends explicitly on the local density of the inhomo- 
geneities: 
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(15) 
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The mixing pressure is negative for both geometries at all 
volume fractions. This means that the metallic phase is 
under "tensile stress" due to the long-range interaction. 
Thus the equilibrium density of the metal is lower than 
the density predicted by MC in the neutral case. For 
a discussion on the stability of a fractionated electronic 
fluid subject to a negative pressure see the Appendix B 
of Ref.lM 

From the definition of the pressure one finds for the 
B-phase chemical potential: 



Ms = 



Jb (n B ) 



Pm 

n B 



The last term is the contribution due to the frustrating 
forces. For the global chemical potential, /i = df /dn, we 
find: 

_ Ib (n B ) m ix 
n B 

with the mixing contribution given by: 

M mi * = ^f ^rPm + ±e m [n B (n),v(n)]. 

We have provided the general equations for a FPS be- 
tween an incompressible phase at ua — and a com- 
pressible phase. To proceed further we need an explicit 
expression for the bulk free energy. 
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FIG. 2: Top: The free energy of the uniform metallic phase 
/h (n) compared with the FPS-state in the case of drops 
[fo (n)] and stripes [fs (n)] for A = 0.02. Bottom: Expanded 
scale in the density range where the geometries compete. 



A. Parabolic approximation for 2D metal free 
energy 

Now we solve our problem expanding the free energy 
of the B-phase in a parabolic approximation around the 
density at which, in absence of the long-range Coulomb 
interaction, the system would experience the transition 
to the phase-separated state. In other words we will label 
n B the density which satisfies the equation p B (n B ) = 0, 
f B the corresponding free energy f B = f B (n B ^j and /x 
the corresponding chemical potential. These quantities 
are related by f B = fJ.on B . 

Using a Taylor expansion we can then write the B- 
phase free energy as follows: 



f B (n B ) = fJ-a n B + 



1 



2k b 



(n B - n° B ) 2 



(16) 



where k b = (d 2 f/d 2 n) 1 is proportional to the com- 
pressibility of the metallic phase. 
The pressure of the B-phase is now: 



PB {n B ) 



2 K B 



(17) 



We introduce a dimensionless parameter A that mea- 
sures the strength of the frustration due to the mixing 
energy effect. A can be defined as the ratio of the charac- 
teristic mixing energy (obtained from Eq. (JSJ without the 



geometric factor u(y)) to the characteristic phase sepa- 
ration energy gain ~ n B 2 / Kb- 



A = 



kb 



(n 



\2T d 



Mb) 



so 



A coincides apart from a numerical factor with the pa- 
rameter introduced in Ref. f° r d = 3. In the follow- 
ing we measure the pressure and free energy densities in 
units of the characteristic phase-separation energy gain 
(n B ) 2 /ks and all the densities in units of the charac- 
teristic MC density n B . Thus we define n' B = riB/n B , 
n' = n/n° B , p' B — pb^bI^b > ^ ne uniform B phase en- 
ergy density fg = fBKBln B2 and the phase separated 
energy density Jd,s = f K B/n B2 . In our parabolic ap- 
proximation the free energy densities read: 

f„{n') = M0 n'^ + ~(n'-lf (18) 
n B z 

Jd,s [n , n B ) = Mo n — + — - [n B - 1) + 

XnBU [y B 

The local density n' B is determined by solving Eq. i|13|) 
with the left and right hand side given by Eq. <|17[) and 
Eqs. (|14fl and (|15fl respectively. 
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FIG. 3: Same as Fig-Ufor A = 0.5. 



To decide the most stable geometry for the FPS and 
its stability against the uniform phase we have to com- 
pare the expressions in Eqs. I|18l) for different strengths 
of the LRC interaction. The first term in the free ener- 
gies represents the MC free energy which is equal for all 
states and can be eliminated. 

In Figs. 12151 we plot the free energies with this term 
subtracted. We show the free energies of the uniform and 
the FPS-state in the presence of drops (/d) and stripes 
(fs) for two different values of the parameter A. For small 
A (Fig. |2J the striped state introduces a lower energetic 
cost with respect to the drop geometry at intermediate 
global densities, while at low and high densities the two 
geometries compete. Increasing the global density one 
finds a first "geometrical" transition from a droplet state 
to a striped one at low density (left bottom panel) and 
a second "geometrical" transition from stripes to drops 
at high density (right bottom panel). For large A the 
situation is different (Fig. |3J.The striped state results 
stable for both intermediate and high density and the two 
geometries compete only at low density (bottom panel) 
where there is a "geometrical" transition from the drops 
to the stripes as in the previous case. 

Note that the free energies in the FPS-state for both 
geometries have a negative curvature and thus exhibit a 
negative electronic compressibility. This, however, does 
not imply necessarily an instability since the background 
contribution to the inverse compressibility, which is very 
large and positive^, should also be taken into account. 



The range of global density where the FPS-state is 
stable depends upon A. Furthermore inside the FPS-state 
one finds that the stable geometric arrangement changes 
both with the global density and the strength of the LRC 
interaction. This leads to construct a phase-diagram in 
the n — A plane (see Fig. HJ. Given our initial choice 
we have two uniform phases, the metallic phase at high 
global densities and the insulating phase at n = 0. For all 
A the two uniform phases are separated by FPS-states. 
The global density range of stability of the FPS, which for 
the unfrustrated case in absence of the LRC interaction 
is determined by < n' < 1 (i.e. < n < n° B ), shrinks 
increasing A. This is clearly due to the tendency of the 
long-range Coulomb interaction to stabilize the uniform 
state as in 3D systems^. Close to the insulator one finds 
metal drops and close to the uniform metal one finds 
circular voids in the metal. We will call the latter "the 
bubble state" . Finally at intermediate densities one finds 
stripes. Increasing A, the bubble state stability range 
shrinks and disappears above a value A* ~ O.f so that 
entering from the metallic uniform phase the FPS-state 
is made of stripe inhomogeneities. On the contrary the 
low density metallic drop state persists at all A. 



1. Properties of the frustrated phase separated state 

We now discuss explicitly the physical properties of 
the FPS-state and the order of the various transitions. 
For the stripe state the relation for the local density as 
a function of global density is obtained in terms of v 
in a parametric form. Then the volume fraction is ob- 
tained from the neutrality constraint. When the inho- 
mogeneous phase has a drop geometrical arrangement 
(metal drop state and bubble state) the local density 
in the metal takes a particularly simple form. Solving 
Eqs. lO ,10,0 with respect to n B we derive the B- 
phase local density in term of the global density: 




(19) 



Obviously in the case A = one recovers the MC solution 
tib = n%. Contrary to the MC the local density of the 
metallic inhomogeneities depends explicitly on and is a 
decreasing function of the global density n of the system. 
Increasing the strength of the long-range interaction A, 
this effect becomes stronger. 

This behavior of the local density versus global den- 
sity can be detected from local probes (like NQR and 
NMR). Physical quantities, which depend on the local 
density, will show an unexpected behavior: they will re- 
spond to an increase in density as if the density were 
decreasing. This effect has been discussed in connection 
with the Curie temperature in manganites^^ 

For the metal drop state at any A, ns approaches the 
MC value n B going towards the homogeneous insulating 
phase, i.e. n — > (Fig. |SJ). Furthermore the transition 
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FIG. 4: Top: The phase-diagram in the n — A plane. Between 
the two uniform phases one finds the FPS-states. The tran- 
sition from the insulator to the metal drop state is second 
order at all A (see Sec. IHI A II . while the transition between 
the metal and the striped FPS-state at A > A* is first order 
as well as the "geometrical" transitions. Bottom: Expanded 
scale for A ~ A*. One finds a second-order transition from 
the metal to the "bubble" state. 



from the FPS to the insulating-phase is second order like, 
since the volume fraction goes continuously to zero with 
sloop 1 [y ~ n') (Fig. EJ). From the same figure we see 
that for A = 0.05 < A*, the system goes from the uniform 
metal to the bubble state with the volume fraction for 
the insulating phase (given by 1 — v) which goes contin- 
uously to zero at the transition. Also the local density of 
the metallic regions is continuous (upper curve in Fig. [SJ) 
indicating that the transition is second-order like. 

For A > A* the region of stability of the bubbles close 
to v ~ 1 disappears and the striped FPS-state appears at 
the transition with insulating inhomogeneities that have 
a finite volume fraction (Fig. |BJ|- From these properties 
one can conclude that the transition is first-order in this 
case. The first order character of the transition can be 
understood from the behavior of the mixing pressure for 
the stripe case. From Eq. (|15|l one notes that in this case 
the mixing pressure has a divergence at = n which 
cannot be reached. Therefore one can not go continu- 
ously from the uniform metal to the stripe state. 
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n/n B 

FIG. 5: The local density behavior vs. the global density for 
A = 0.05, 0.5. In correspondence to the "geometrical" tran- 
sitions and at the transition from the metal to the striped 
state tib shows a discontinuity reflecting their first-order na- 
ture. The segments close to the origin corresponds to metallic 
drops. The small segment in the upper curve close to the full 
line is the bubble state. The rest of the dashed curves corre- 
spond to stripes. The full line is the uniform metal. 



Inside the FPS-state one can see that the high density 
"geometrical" transition from the "bubble" state to the 
striped one for A < A* and the low density "geometrical" 
transition from the latter state to the metal drop state 
have a discontinuity of both the volume fraction and the 
B local density [Figs. reflecting the first-order na- 

ture of the transitions. This is expected since at least in 
our approach one can not continuously deform drops to 
get stripes. 

For all the FPS-states, increasing the strength of 
the long-range interaction the metallic density decreases 
(Fig. [SJ in order to minimize the mixing energy of the 
inhomogeneous state which is oc Xns- At the same time, 
increasing A (Fig.[SJ), the volume fraction has a growing 
rate larger than in the ordinary MC case. 

The order of the transitions can be also checked if one 
looks to the chemical potential /i = df/dn (Fig. 0. In 
the range of stability of the FPS "bubble" state (small A) 
the chemical potential at the critical density has a cusp 
indicating the second-order nature of the transition. In- 
creasing A to values greater than A* the "bubble" state 
disappears and the chemical potential has a discontinuity 
at the transition to the striped FPS-state. A similar dis- 
continuity in \x is also obtained in correspondence to the 
"geometrical" transitions consistent with the fact that 
these transitions are first-order like. Notice that for the 
insulator, which is incompressible, the chemical potential 
is not defined. 

Finally we analyze the size Rd of the inhomogeneities 
and the size R c of the cells in which the system is di- 
vided. One can introduce a length scale Id which fixes 
the natural units for the characteristic size of the cell. 
This scale correspond to R c (Eqs. iJB}, Q) evaluated at 
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FIG. 6: Top: The volume fraction of the B-metallic phase in 
the FPS-state vs n at different values of A compared with the 
ordinary Maxwell construction analysis. The full horizontal 
line corresponds to the uniform metal. For A = 0.05 the 
small segment close to the uniform metal is the bubble state. 
Close to the origin for both values of A one has the drop 
state. The rest is in the stripe state. Bottom: Expanded 
scale at low density to show the abrupt change of v due to 
the "geometrical" transition from the metal drop state to the 
striped state as density increases. 
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FIG. 7: The chemical potential behavior in the FPS-state. 
When the "bubble" state exists, at the transition from the 
metal to the FPS-state n exhibits a cusp point (2 nd order 
transition) while for A > A* the transition to the striped state 
implies a chemical potential discontinuity (1 st order transi- 
tion) 
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FIG. 8: The cell radius and the inhomogeneity radius vs the 
global density for small A. The divergence in R c are in corre- 
spondence to the transitions from the FPS-state to both the 
metal and the insulator. 



the MC-density n B dropping all the geometric factors: 

(20) 
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Another important length scale is the screening length: 



L 



2ire 2 kb 



(21) 



Inserting the last equation in the definition of A one finds 
that A = ld/(2nl a ). Similar length scales where defined 
in the 3D casej^i in particular we can write the screen- 
ing length for general dimension (neglecting factors of tt) 
l s ~ (eo/(e 2 KB)) 1 ^ d ~ 1 - In 3D systems Is gives an upper 
limit for the size of inhomogeneitiesji^ii&2i The nature 
of screening is quite different in 2D systems^ As it will 
be clear in Sec. IIVI l s still plays a fundamental but does 
not limit the size of inhomogeneities. 
With the above definitions one has: 



R, 



= 2tt\I s — 



u{v) 



Stripes 



n' B u (u) 



Drops 



(22) 



(23) 



For the stripes the half-width of the metallic regions is 
given by Rd = R c v. For the drop geometry it is more 
convenient to define a symmetrized radius Rd = R c [v(l — 
v)] 3 that smoothly interpolates between the radius of 
metallic drops at small v and the radius of bubbles for 
v close to 1. The above assumptions lead to the simple 
expression, 



Rd = V3~7T 



AL 



V3 



(24) 
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For small n' or small A one has an almost constant be- 
havior R d ~ 0.9l d [c.f. Figs. 1519]. 

For A < A* bubbles appear in the metal with a diver- 
gence of the cell radius while the drop radius remains 
finite (see Fig. |SJ. That is at threshold bubbles appear 
suddenly with a finite size but the transition is second 
order because they are infinitely far apart. This is remi- 
niscent of the transition in a type 11 superconductor as a 
function of field at H c 1 , which according to GL theory is 
second-order although normal state "drops" (the vortex 
core) have a finite radius £^1 

For A < A*, -Rd is of the order of l s . If the latter is 
of the order of the interparticle distance, a mesoscopic 
treatment may be problematic. On the other hand in- 
creasing A, Rd grows to values that are much larger than 
Is (Fig®. This represents the main difference respect to 
3D systems for which it was demonstrated that l s is an 
upper bound for Rd at any A. In the next section we will 
show that this difference stands on the different role of 
the LRC interaction and screening in 2D and 3D systems. 

The divergence of the cell radius disappear when A is 
greater than A*, since at the first-order transition to the 
striped state the cell radius stays finite. 

Increasing A at fixed global density one finds that R c /l s 
increases. This behavior is easy to rationalize if one con- 
siders an increase of the surface energy in such a way 
that A increases while keeping l s constant. In this case 
the system prefers to make domains with longer period- 
icity to reduce the surface energy. From Fig. we see 
that metallic stripes become narrower as the insulator is 
approached. 



IV. LOCAL DENSITY APPROXIMATION 
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FIG. 9: Top: R c and R d vs. n for A=0.08. The FPS-state 
appear with a finite cell radius. The typical size of the inho- 
mogeneities is of the order of Is- Bottom: Same of top for 
A = 0.48. The sizes of the domains and the cells are much 
greater than Is. 



In this section we want to relax the uniform density 
approximation. For simplicity as in Sec. IIIII we restrict 
ourselves to study a phase separation between a metallic 
phase B and an insulating phase A with energy density 
1a = o. 

In order to account for the spatial dependence of the 
local density, we take the total free energy as a local den- 
sity functional. Indicating with Vb the B-phase domain 
volume and with V the total system volume the total free 
energy can be put as: 



F = / f B (n B (r))d 2 r + a^ AB + ^- 

JreB Ze 



[n B (n) - n] [n B (r 2 ) - n] d 2 r 1 d 2 r2 



\n - r 2 \ 
(25) 



where n is the global density, while T, AB is the total 
interface surface. We still assume a sharp interface due to 
short range forces with the parameter a parameterizing 
the surface energy. Eventually the surface energy term 
could be replaced by a gradient term to consider soft 
interfaces. 

The constraint of charge neutrality of the system reads: 



/ 



tib (r) d r = nV 



(26) 



We minimize the total free energy functional with respect 
to the local density and include the above neutrality con- 
straint via a Lagrange multiplier \x e . We then obtain an 
equation which tell us that the electrochemical potential 
is constant: 

Ms(r)-— <t>{r)=iie VreB (27) 

where (f> (r) indicates the electrostatic potential gener- 
ated by the charge distribution [n B (r) — n]. This equa- 
tion has to be solved together with the neutrality condi- 
tion Eq. dp . 

Assuming again for the B-phase the parabolic free en- 
ergy [Eq. l(Po|) ]. one obtains an equation relating the local 
density to the potential in terms of the compressibility: 



n B (r) -n° B = k b 



—cj) (r) - (no - p, e ) 
£o 



VreB (28) 



10 




X' 



FIG. 10: The excess charge density profile n^ =0 0') In - 1 
in the direction perpendicular to the stripe for v = 1 /4. x' is 
the component of r / R d 



or in terms of the screening length 



(r)- — (/Xo-Me) 
e 



27reZ s [nzj (r) - n%] VreB (29) 



In the limit of infinite compressibility i.e. zero screen- 
ing length, the electrostatic potential is constant on the 
metallic regions and therefore tib (t) correspond to the 
distribution of a metal for which the 3D Laplace equa- 
tion V <fi — is supplemented by the boundary condition 
4> = const on the domains. Here we are assuming that 
h/ Rd — ► but we are not making assumptions about the 
other parameters so that A ~ ld/h may remain finite. 

In this limit the problem can be solved analytically for 
the stripes geometry. In fact, in this case the Coulomb 
potential is calculated by using the Schwarz-Christoffel 
conformal transformations^ The B-phase density spa- 
tial dependence comes out to be: 



cos ■ 



(30) 



sin 



where x' indicates the a; component of the dimensionless 
coordinate r' defined by r' = r/Rd- In Fig. 1)10(1 we 
show the spatial dependence of the excess charge density 
n l g~° (a/) /n — 1 on one unit cell. 

Keeping in mind that the charge in the stripe region 
is undercompensated by the background, one finds that 
at this level of approximation the metallic stripe behaves 
as a macroscopic charged metallic strip. The density ac- 
cumulates on the border of the stripe and decays as a 
power law towards the center. 

It is interesting to compare this configuration with the 
three dimensional analogue, that is metallic layers locally 
undercompensated by a uniform background. In that 
case the charge is localized at the surface of the layer 
and decays as exp(— x/l s ). Although in both cases l s is 
the length scale below which macroscopic electrostatic, 



i.e <j) = const., is not any more valid, the solutions for 
x » l s are dramatically different in the two cases. For 
3D metals the local density far from the surface is forced 
to be equal to the background density which in turn is 
equal to the global electronic density i.e. ns — On 
the contrary on the present 2D case the local density far 
from the surface can be quite different from the global 
density. 

In Ref. I3l3l it was shown that for small Is /Id in the 
3D case mesoscopic inhomogeneous phase separation is 
forbidden. In fact the system gains phase-separation en- 
ergy when the local density differs from the global density 
value. Instead electrostatic forces ub = n, except for the 
microscopic length Is, and the system remains uniform. 
In the 2D case the inhomogeneity is able to gain phase 
separation energy in all the region where ns ^ n which 
is not limited by a microscopic length. 

A related issue in 3D is that in general inhomogeneities 
can not have all linear dimensions larger than l s . In the 
2D case this does not impose any constraint because the 
linear size perpendicular to the layer is from the outset 
smaller than l\ D . Therefore inhomogeneities can have an 
unbound size in the plane. 

Notice that this different behavior between 2D and 3D 
systems is due to the difference in the charge profile far 
from the surface. Of course the divergence of the elec- 
tronic density at the surface of the stripe is unphysical 
and will be cutoff by a microscopic length (see below) but 
this does not affect the behavior of the charge far from 
the surface that is essential for our argument. 

In the related problem of a two-dimensional electron 
(hole) gas in MOSFET devices the background is not 
rigid but is provided by a metallic gate in the form of 
mobile holes (electrons). In this case the background 
can relax and other considerations enter into play^iS 
In the case in which the separation between the gate and 
the e-gas is much larger than R c and Rd, the electric 
field produced by the electron gas inhomogeneities will 
be perpendicular to the gate plane and this will prevent 
the background to relax. This is illustrated in Fig. ITU 
where we show the equipotential lines in the plane per- 
pendicular to the stripe. We see that at a distance ~ R c 
from the stripe the equipotential lines become parallel to 
the stripe surface. In this case the background behaves 
as rigid and our results apply. In other words the back- 
ground can follow the inhomogeneities of PS only when 
the inhomogeneities are comparable in size to the sep- 
aration between the e-gas and the gate. The crossover 
between these two possible scenarios will be considered 
elsewhere 

To further clarify the role of Is, now we calculate the 
charge density profile for Is ^ adapting the method 
of Ref. 0* The unphysical divergence of the metallic 
density will be removed. We restrict to the stripe geom- 
etry as above but the same method can be used for other 
geometries. Consider the 3D Poisson equation 



(r,z) = -47175 (r,z) 



(31) 
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FIG. 11: The equipotential lines in the plane perpendicular 
to the stripes for the same parameters as Fig. 1101 



with z the distance from the plane and the 3D charge 
density denned as p (r, z) — p(r)S (z) with p(r) the 2D 
in plane density. To solve Ea. (|31[l . we perform the Fourier 
transform in the x direction using the fact that the solu- 
tion is periodic in 2R C : 



:3D 



<f> 3D (q,z) 



Rc 



dxe- iqx 4? D (x, z) 



(32) 



(and similarly for the charge distribution). Since there is 
no dependence in the direction of the stripe we dropped 
the y coordinate. The Poisson equation is then: 



d 2 r D (q,z) 
dz 2 



q 2 ^ D (q,z) = -47Tp(q)S(z) (33) 



with p(q) defined as the Fourier transform of the in plane 
density 

p(r) = -e [n B (r) - n] . 

For z the solution for the potential is 4> 3D (q, z) = 
3I) (g,O)e-l9ll 2 l. Integrating Eq. (33) in a small interval 
one obtains the boundary condition: 



\q\<t>{q) = 2irp(q) 



(34) 



with the Fourier transform of the in-plane electrostatic 
Coulomb potential defined as cf>(q) = <fi 3D (q, 0). The 3D 
boundary condition Eq. I|34|l looks formally like an effec- 
tive 2D Poisson equation. 

Eq. (|28|l or Eq. I|29|l determines how the charge re- 
sponds to the potential and Eq. Q34JI determines how 
the potential is generated by the charges. Both equa- 
tions must be solved self-consistently in order to find the 
charge distribution. 

Using the superposition principle both charge and po- 
tential can be written as the sum of the terms evaluated 



above for infinite compressibility (i.e. for l s 
correction, which we wish to compute: 



<j) ls=Q (r) + Scp (r) 

' B 



tj>{r) 

n B (r) = n 1 ^ (r) + 6n B (r) 



0) plus a 

(35) 
(36) 



The correction S(f> (r) satisfy the effective Poisson equa- 
tion: 



\q\6(f>(q) = -2Tre6n B (q) 



(37) 



The unknown Lagrange multiplier p, e has to be deter- 
mined by fulfilling the neutrality condition and can also 
change as l s is increased from zero: 



Me = m; 
Eq. (|2*5|> can be put as 

5n B (r) = 1 5(j>(r) 



5n e 



2ne Is 
Vr e B 



-°(r) 



(38) 



where we have absorbed the Lagrange parameter in the 



constant n e = n° B 
8n B (q) 



2Tve 2 l a 
1 



In Fourier space we get: 
sin [(q - q') R d ] 



2ire IsRc 



(q - <i') 



Js=0 



(?) 



2sin(qR d ) 



(39) 



Eq. I|37|l and Ea. H39l) are a closed system since the quanti- 
ties with Is — are known from the previous treatment. 

In the case R c ,Rd >> Is one can substitute J2 q ~ * 
2R C J dq/{2it) and make the approximation: 



2 sin [{q-q')R d 



(q - 9') 

Using Eq. (|37|l we obtain: 

2n e sin(qR d )/q 



2n5(q-q'). 



Sn B (q) = 



Js=0 
B 



l + Ws\q\) 



(40) 



(41) 



In the limit q — * we obtain Sn B — > 0. The uniform 
component of the charge does not change and therefore 
<5/i e = and n e = . We have evaluated the above 
expression via a discrete Fourier transform in the limit in 
which the stripes are far apart (R c /R d — > oo i.e v — > 0). 
This correspond to solve the problem for a single stripe. 
The electronic density at Is = can be put as: 



=° [x>) 



1 



™ y/1 - (a/) 



IY2 



In Fig. E| we show the total electronic charge density 
for ls/Rd = 0.03. The main difference with respect to 
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FIG. 12: Top panel: Comparison between the electronic 
charge density at ls/Rd = and ls/Rd = 0.03 evaluated 
at v = 1 and n = 1. Bottom panehExpansion near the inho- 
mogeneity surface 



the Is = case is that the unphysical divergence of the 
density at the stripe surface is removed and the density 
tends to a finite value at the stripe boundary. Notice 
that at distances from the interface larger than lg the 
charge density at Is ^ practically coincides with the 
one at Is = as expected. Our previous conclusion re- 
garding the absence of upper bounds for the size of the 
domains remains unchanged. In the entire domain the lo- 
cal density differs from the average density and therefore 
the phase separation energy gain comes from the whole 
domain and not from the electric field penetration depth 
as it happened in the 3D case. 

Finally one has that the electronic potential is given 
by: 



2?re 



1 



2 sm(qRd) 



(42) 



The above electrostatic potential includes the two- 
dimensional screening with lg as a screening parameter. 
It is of the same form of the screened potential discussed 
in Ref.lll 



V. DISCUSSION AND CONCLUSION 

In this work we have considered the problem of meso- 
scopic FPS in 2D electronic systems frustrated by the 
LRC interaction and the surface energy. In particular 
we concentrated on the problem of coexistence between 
a 2D metal and an insulator. 

We have supposed that the inhomogeneous state is re- 
alized with two different geometric arrangements: disks 
of one phase into the other and a state with alternat- 
ing stripes of metal and insulator. The first arrangement 
comes in two different flavors: bubbles of insulator hosted 
by the metal for v ~ 1 and metal drops hosted by the 
insulator for v ~ 0. 

We have defined a parameter A which specifies the 
strength of the frustration and is given by the ratio of the 
characteristic mixing energy to the characteristic phase 
separation energy gain in the absence of frustration. 

As in the 3D case, we have found that frustration tends 
to extend the global density region where the uniform 
metallic phase is stable. Within our approximation, be- 
low some critical value of A, we have found that as the 
global density is reduced there is a second order transition 
from the homogeneous metallic state to an inhomoge- 
neous state with insulating bubbles in the metal. Above 
this value of A the bubble state is never stable and the in- 
homogeneous transition leads to the striped state with a 
first-order transition. Inside the inhomogeneous stability 
region by decreasing the global density first-order "topo- 
logical" transitions induces changes in the geometry of 
the domains. 

In Sec. IIII A II we have analyzed the inhomogeneous 
state properties and we have found that the size of in- 
homogeneities is not forced to have one linear dimension 
smaller than the screening length l s in sharp contrast 
with the 3D case. This difference stands on the qualita- 
tive different behavior of the screening in a 3D and in a 
2D system. In the first case the screening decays expo- 
nentially from the interfaces whereas in the second case 
it decays as a power law. In this last case it is possible 
to gain PS energy from regions far from the skin depth. 

Inhomogeneities seem to be related to a number of in- 
teresting phenomena like colossal magnetoresistent. Our 
study suggests that inhomogeneities may be favored by 
engineering materials with enhanced anisotropy. For ex- 
ample adding insulating layers in between metallic layers 
should favor mesoscopic phase separation. 

A realistic model of the 2DEG would require consid- 
ering a metallic phase with negative compressibility as 
observed in the uniform phase. The effects of disorder 
should also be included4i Therefore here we limit our dis- 
cussion to some generic qualitative features. The stabi- 
lization effect of the long-range Coulomb interaction has 
been observed on the 2DEG where the state with nega- 
tive short-range electronic compressibility has been show 
to be stable in a certain range of densities ^i^^/S^L^ 
Hani and collaboratoraSiiS observed that the system be- 
comes inhomogeneous at a mesoscopic scale when the 
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system is at the verge of the metal-insulator transition 
suggesting a connection between the two phenomena. 

Another interesting finding of Hani and 
collaboratorsSiiS is that coming from the metallic 
side the compressibility has negative spikes close to 
the transition. When integrated those spikes imply a 
decreasing step in the chemical potential as a function 
of density. This behavior is similar to what we find 
for the transition from the inhomogeneous state to the 
metallic state as the density is increased in the case in 
which the transition is first order (c.f. Fig.[7J). The main 
difference between our result and the experiment is that 
in our case there is a single large step at the transition 
whereas Hani et al. find many small steps around the 
critical density. This behavior is easy to rationalize if 
one considers that in the presence of disorder the density 
will not be uniform in all the sample. A distribution 
of critical densities and fragmentation of the large step 
in many small steps will be naturally produced. The 
minimum step size will be given by the appearance of a 
single drop. 

We have restricted ourselves to one isolated plane but 
we expect the same physics to be true for weakly coupled 
planes as found in some cuprates, nickelates manganates 
and other strongly correlated materials. As mentioned 
in the introduction the situation in manganites is quite 
complex due to the large variety of competing phases. In 
layered cuprates one finds striped charge and spin den- 
sity waves at small doping and a Fermi liquid at large 
doping. From the behavior of the mean-field energy^ it 
is quite natural to expect phase separation among the 
stripe state and the overdoped Fermi liquid. To make 
a rough estimate of the drop size we identify the stripe 
state as the pseudogap phase and treat it as incompress- 
ible. If we assume that the surface energy is of order 
a ~ J/a, with J the superexchange constant and a the 
lattice constant, the drop radius can be written as [c.f. 
Eq. 

a_( e aJ \ 
Kd d Ax\2a Ry) 

where ao is the Borh radius, Ax ~ 0.1 is the range of 
doping where one expects phase separation according to 
MC. Using J — Q.OlRy, a — 7.2ao, eo = 5, one gets 
2Rd ~ 8a ~ 30A which compares well with the inho- 
mogeneities found in Refll3l. This experiments reveals 
superconducting like inhomogeneities of roughly 50A di- 
ameter (playing the role of the compressible phase) em- 
bedded in insulating like regions. A similar rough esti- 
mate gives A ~ 1. 

The constraint on the maximum size of inhomo- 
geneities in 3D systems makes the conditions for stability 
of a phase separated state very stringent i 3 ^ In 2D system 
instead we have found that this constraint does not apply 
and frustrated phase separation at the mesoscopic scale 
is much more favorable. This may be the reason why in- 
homogeneities are often found in quasi two-dimensional 
electronic systems. 
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APPENDIX A: ELECTROSTATIC ENERGY IN 
THE UNIFORM DENSITY APPROXIMATION 

Here we compute the explicit expressions of the elec- 
trostatic energy within the UDA for the two geometries 
considered. We start by considering regular arrays of in- 
homogeneities and dividing the system in Wigner-Seitz 
cells in such a way that each cell is globally neutral. The 
electrostatic energy is iV<j times the Coulomb energy of 
one cell. 

As in Ref . H5 we take the density profile inside each cell 
as —e{riB — tia) inside the B-phase inhomogeneity of 
radius Rd compensated by the background charge density 
e(n- n A ). 

In 3D the computation of the Coulomb energy is facil- 
itated by the use of Gauss theorem^ In the present 2D 
case with the full 3D Coulomb interaction Gauss theorem 
is not useful so we use the explicit expression: 53 

e e i = - — tjt I d 2 r (p d + p b ) (fa (r) + c/) b (r)) (Al) 
2 £o V d J Vd 

where Vd is the 2D volume of one cell, pd = 
— e {jib — tia) is the inhomogeneity charge density and 
Pb = e (n — tia) is the effective background charge den- 
sity. We need to evaluate the Coulomb potential <pd 
generated by the inhomogeneities of all cells and the 
Coulomb potential 4> b generated by the whole background 
charge density. The latter corresponds to the Coulomb 
potential of an infinite uniformly charged plane and can 
be taken as constant in the plane. Because of the global 
neutrality the value of the constant does not affect the 
electrostatic energy and we can take 4> b = 0. 

The inhomogeneity Coulomb potential is the sum of 
the Nd contributions from each cell. We separate the 
contribution of the cell where we are integrating from 
the contribution of the other cells. In this way Eq. (|A1() 
separates in a self-energy contribution (e^) and an inter- 
action contribution e™*: 

£l = TT \ d 2 r(p d + p b )^(r) (A2) 
2e J Vd 

elf = T- / d2r (P d + Pb )t mt (r) (A3) 
2e Jv d 

1. Stripe geometry 

In this geometry we assume that the system is divided 
in cells of width 2 R c and length L with periodic bound- 
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ary conditions. Within each cell the width of the inho- 
mogeneity with charge density pd is equal to 2 Rd- 

To compute the electrostatic energy we have evaluated 
the two expressions Eqs. I|A2J| . I)A3|1 . The interaction 
potential (f> has been numerically computed truncating 
the sum of the contribution from each cell to a finite 
number of cells: N c . For N c not too large the Coulomb 
energy is asymmetric. However for N c of the order of 
10 2 the Coulomb potential becomes symmetric respect 
to the " phase-exchange" symmetry (A — * B , v — *■ 1 — v) 
indicating the achieved convergence. 

In the limit v — > 0, Rd << R c , <\> mt can be neglected so 
that the total electrostatic energy is well approximated 
in this limit by: 



lilr i eel = — (n>B 



n A f R c 2v 2 [- hxv] 



(A4) 



We can obtain a similar approximation in the opposite 
limit v —> 1 by imposing the "phase-exchange" symme- 
try. This leads to the approximate expression Eq. (J5J 
which interpolates between the v — ► 0, 1 limits. The 
comparison between the approximate and the numeri- 
cal result is shown in Fig. ^ in which are reported the 
U functions considering the different electrostatic ener- 
gies. We see that the approximation to keep only the self 
energy term is indeed very good. 



2. Drop geometry 

From the computation in the stripe geometry it is clear 
that in the limit of v — > and v — > 1 the dominant 
electrostatic term is the self-energy. Since the drops are 
zero dimensional objects we expect that the effect of the 
intercell terms to be even smaller than for the case of 
stripes. To check this hypothesis we first estimate the 
interaction energy between the cells in this way: if we 
look to the total system of Nd cells in the limit of small 
volume fraction, the drops can be considered as negative 
point-charges at distance R c since v — > is equivalent to 
Rd << R c and we assume that the drops arrange in a 
Wigner crystal. The electrostatic energy is given byi^ 



Eft* oc -a 



[eR 2 d (us - n A )\ 
R, 



(A5) 



Since v = R 2 ,/ R 2 for the drop geometry and referring to 



the energy per unit volume one has 

s 2 
to 



e eV « — {n B - n A f R c av 2 



The Coulomb self-energy of a cell of radius Rd can be 
evaluated noting that the cell charge density in units of 
(— e) can be written in the Fourier space as^: 



Sn (k) = 2tt (ub — n A ) Rd 



Jx (kR d ) Jx (kR c 
kRd kR c 



where Sn (k) = [rig (k) — n A (k)] 
ergy per unit volume is: 



The electrostatic en- 



^- 2 {n B {k)-n A{ k)) 2 ^ 



61 2£ 7Ti? 2 J ( 27 r) 

Computing the coulombic self-energy per unit volume in 
the Fourier-space one obtains: 



= — (nB -n A ) R c - v 2 + o [y ) 



so 



(A6) 



were we have kept the dominant contribution when v — > 
0. Also for drops one obtains that in this limit the self- 
energy term dominates the electrostatic energy. 

The limit v — ► 1 can be obtained by replacing v — > 
1 — v. This leads to the approximate electrostatic en- 
ergy expression Eq. which as for stripes interpolates 
between the two limits. 

Notice that since the inter drop interaction is negligible 
our computation is independent of the lattice structure 
of the crystal and is also valid for an amorphous config- 
urations of drops. 



APPENDIX B: LIMIT OF SMALL SCREENING 
LENGTH IN THE UDA AND COMPARISON 
BETWEEN THE UDA AND THE LDA 

To compare the LDA and the UDA we use the particu- 
larly simple limit Is — > 0. Formally this can be achieved 
by making Kg -> oo so that the B phase bulk energy be- 
comes a linear function of density. If n° B is kept constant 
the line intersects the energy of the A-phase in this limit. 

We will instead consider the energy density difference 
AE g between the B and the A phase at zero density 
fixed. In this case the MC density n° B = [AE^kb] 1 / 2 
diverges. 

To make things less abstract we can consider the fol- 
lowing example: classical electrons with a short range 
attractive interaction in the lattice and at low tempera- 
ture. The attraction stabilizes a crystal at high density. 
The energy per unit volume of the uniform crystal is 
= —zvno/2 where z is the coordination number v is 
the short range attraction and no is the density of the 
incompressible crystal phase. 

For small global density the attraction can be ne- 
glected. Electrons form a uniform "metallic" phase. 
Since the electrons are classical (no tight binding hop- 
ping term) the chemical potential becomes independent 
of the density (thus the compressibility is infinite) and 
taken to be zero. 

In the absence of long-range interaction and for elec- 
tronic densities n e < no this system phase separate into 
the high density crystal phase and the zero density empty 
phase. In the presence of the long-range Coulomb inter- 
action separation between the high-density crystal phase 
and a low-density metallic phase becomes possible. 
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In order to keep notations consistent with the previous 
section we consider the hole density n — n Q — n e . The 
incompressible A phase is then at n = and has energy 
f A < and the infinitely compressible metallic phase has 
n > and energy = 0. 

The free energies of the homogeneous state and the 
PS-states (drop and stripes geometries) can be put as: 



For the stripes geometry the energy density reads: 



f H (n) 
Jd,s 0, v) 







(i + 



£0 



(Bl) 

n^l (B2) 



Minimizing this expression respect to v one obtains the 
optimum v value. 

Now in order to test the accuracy of the UDA we will 
compare this approximation with the LDA derived in Sec. 
IIVI For simplicity we restrict to the case of separation 
between and incompressible phase and an infinitely com- 
pressible phase. 

The energy can be evaluated in the LDA using the spa- 
tial dependence of the electron density found in Sec. IIVI 
for Z s = [c.f. Eq. ■ 

As discussed above we assume that Ja is a negative 
constant and 

Sb K (r)] = 

With these conventions the LDA free energy functional 
reads: 

F = V A f A + a^AB + ^- I n4> ls=0 (r) d 2 r (B3) 
2eo Jv A 

where we have used that in the metallic phase domains 
the electrostatic potential for Is — is constant and thus 
give no contribution to the Coulomb energy. 



f s = (l-v)f A+ n 2 R, 
4e 



V 2 Rc 



Here ulda is given in terms of the potential at Is = 0: 

6' s=0 (x) dx 



ULDA W) 



,2 ^ n 



e nR c 



where x is the dimensionless coordinate x/R c and the 
potential at Is = reads: 



sin 



= 4i? c en cosh- 



Minimizing the mixing energy e m 
the cell radius R c one finally finds: 



(B4) 



e e i + e a respect to 



.ULDA {v) 



f S {%v) = {l-v)f A + 

V £ o J v 
The difference of energy between the PS and the homoge- 
neous states can be put in the same form as in the UDA 
Eq. I|B2(1 . The only difference with the UDA is encoded 
in the function ulda- In Fig. (TJ we compare the two ge- 
ometrical u function parameterizing the mixing energies. 
The LDA function implies a lower mixing energy since 
we are relaxing the uniform density constraint. This re- 
laxation energy gain however is small demonstrating the 
accuracy of the UDA approximation for thermodynamic 
quantities as it was found also in the 3D case^ 
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